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We present a full phase diagram for the one-dimensional (ID) to three-dimensional (3D) crossover of the 
Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) state in an attractive Hubbard model of 3D-coupled chains in a har- 
monic trap. We employ real-space dynamical mean-field theory which describes full local quantum fluctuations 
beyond the usual mean-field and local density approximation. We find strong dimensionality effects on the shell 
structure undergoing a crossover between distinctive quasi- ID and quasi-3D regimes. We predict an optimal 
regime for the FFLO state that is considerably extended to intermediate interchain couplings and polarizations, 
directly realizable with ultracold atomic gases. We find that the ID-like FFLO feature is vulnerable to thermal 
fluctuations, while the FFLO state of mixed 1D-3D character can be stabilized at a higher temperature. 



The interplay between fermion pairing and magnetism is 
at the heart of understanding strongly correlated systems 
ranging from unconventional superconductors and ultracold 
gases to neutron stars and quarks.' BCS-type superconduc- 
tivity is suppressed by a large magnetic field exceeding the 
Chandrasekhar-Clogston limit. However, it has been pro- 
posed as a paradigm of superconductivity in high magnetic 
fields that it is possible for superconductivity and magnetism 
to coexist with exotic pairing mechanisms.^ ^ The Fulde- 
Ferrell-Larkin-Ovchinnikov (FFLO) state would arise with 
this interplay,^ ^ but it still remains elusive in spite of indirect 
experimental evidence observed.^ The FFLO state is charac- 
terized by the Cooper pair carrying finite momentum causing 
a spatially modulated order parameter. One peculiar feature 
of this exotic phase is that apparently its stability is largely 
affected by the dimensionality of the system. It turns out that 
the three-dimensional (3D)-FFLO state occupies a thin area 
of the mean-field phase diagram^ though the signature can be 
stronger in the systems that support nesting such as optical 
lattices^ ^ and elongated traps.™ Indeed, only phase separation 
was observed for ultracold gases in 3D traps. ' ' On the other 
hand, in the exact solution of a one-dimensional (ID) system, 
the FFLO character appears at any finite spin-polarization,-^ 
while long-range order cannot exist in ID. The experimeniSl 
in ID was consistent with the FFLO theory, although the state 
remains unidentified. 

A natural question arising is whether one can combine the 
promising ID-FFLO features and long-range order provided 
by higher dimensions. This has been considered in previous 
mean-field'^ and effective field theorjl^l studies for coupled 
continuum- ID gases, and for the Hubbard ladder.'^ The trap- 
ping potential is essential in ultracold gas experiments, and 
thus in spin-polarized systems, one can expect a shell struc- 
ture of different phases along the trap. Therefore, beyond the 
local density approximation, the inhomogeneous superfluid 
and normal phases need to be treated in a unified framework 
by including full local quantum fluctuations. Neglecting lo- 
cal quantum fluctuations creates an apparent bias in favor of 
the superfluid state and against the normal state.'*' Here, us- 
ing a real-space dynamical mean-field theory (DMFT), we in- 
vestigate the 1D-3D crossover problem within the anisotropic 
Hubbard model in a trap. 

The dimensionality effect to the FFLO state that we are 



considering here differs from the two-dimensional (2D)-3D 
crossover studied in layered superconductors.'^ There, the 
quasi-2D character minimizes the orbital pair breaking effects 
in a magnetic field, and the Zeeman effect may then lead to 
a d-wave FFLO state. Here we consider s-wave pairing, and 
orbital effects are absent. Note that a particle-hole transfor- 
mation can bring in another interesting perspective by map- 
ping the FFLO state to the striped phase of the doped repul- 
sive Hubbard model.-" This emphasizes the importance of the 
FFLO state in the general context of high-Tc superconductiv- 
ity. 

We perform a real-space variant of DMFT calculations on 
the attractive Hubbard model of 3D-coupled chains, 

ila {II') i<y 

il ila 

where cli^icua) creates (annihilates) a fermion with spin 
a =t,i at site i of chain I, the density fiu^ = cli^cu^r- 
We define the superfluid order parameter as A = — (c|c|). 
Throughout the calculations, the hopping i|| is set to unity. 
The dimensionality is thus tuned by the interchain coupling 
tj_ and varied from ID (tj_ = 0) to 3D (tj_ = 1). The 
chemical potentials fi^ and control the polarization P ~ 
{N^ ~ Ni)/{N^ + ^l) by keeping the total particle number 
+ Ni ^ 120. The harmonic potential trapping the gases 
in a longitudinal chain is given as V^^ = 5 x 10~^{i — 1/2)^. 
The on-site interaction U is selected to be the value corre- 
sponding to unitarity where the 3D two-body scattering length 
diverges.^' The value of U varies from —2.038 {t± = 0.1) to 
— 7.915 (ij^ = 1), which is comparable to the half bandwidths. 

In order to treat an inhomogeneous phase along the 
chain, a site-d ependent self-energy is considered within 
DMFTP^I ^^ I ^^ I With translational invariance in transverse di- 
rections, the on-chain Green's function is written for sites i, j 
and transverse momentum fcj^ as 

[G"^(fc_L;ia;„)]y = [iuJ„(To - ek^as - Si(ia;„)](5,y - /if^-, 

where a;„ denotes the Matsubara frequency, cr is the Pauli ma- 
trix, /I'l is the non-interacting part of the chain Hamiltonian, 
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the 2D energy dispersion e^^ = — 2<j^(cos + cos ky) with 
fcj^ = (kxjky), and the operators are in the Nambu basis. 
In the self-consistency of DMFT, the impurity Green's func- 
tion is obtained as Qi{iujn) = J2k± Gu{k±-,iujn)- While we 
have a local but site-dependent self-energy term, the proce- 
dures can be compared to the chain- dmftP The exact di- 
agonalization method is employed to solve the site-dependent 
impurity problem.'" 

With the on-site interaction U = —4, we have found 
that our calculations qualitatively reproduce the two impor- 
tant characteristics of strongly interacting ID Fermi gases.l^^l^ 
First, at small polarizations, the trap edges are fully paired 
while the trap center is partially polarized. The entire area be- 
comes polarized as the total polarization increases. Second, at 
finite polarizations, the polarized trap center is associated with 
the FFLO state exhibiting a spatially oscillating order param- 
eter 

Clear features of getting away from the ID limit are ob- 
served. For the studied finite interchain couplings, the pairing 
at the edges, namely the first ID feature listed above, is eas- 
ily broken at small polarizations. This indicates that imme- 
diately away from ID, one can observe fully polarized edges 
much before the whole area gets polarized. In contrast, it turns 
out that the other part of the ID features, the emergence of a 
partially polarized center in fully paired vicinities, survives at 
small interchain couplings. Away from but close to the ID 
limit, we typically find a shell structure of polarized edges, 
fully paired shoulders, and a partially polarized center 

In Fig.[T| we present a real-space DMFT phase diagram of 
the Hubbard model of 3D-coupled chains at zero temperature. 
Considering only the phase at the trap center, we find that the 
emergence of the FFLO-type oscillating order parameter oc- 
cupies a large area of polarizations and interchain couplings. 
This wide coexistence area of a finite density difference and a 
finite order parameter extends all the way to the 3D limit, as 
opposed to the mean-field phase diagram on a 3D continuum 
where the FFLO phase occupies only a tiny area.l^ 

Our zero-temperature phase diagram characterizes the 
emergence of three types of shell structures, as shown in 
Fig. [TJb). First, area I shows quasi-3D features in the shell 
structure where the fully paired superfluid (SF) core exists. 
There are FFLO-type oscillations found between the SF core 
and the fully polarized edges. Second, area II shows an in- 
verted sequence: an FFLO core surrounded by fully paired 
shoulders. The edges are polarized in this area, and small or- 
der parameter oscillations also exist at the interfaces between 
the fully paired shoulders and the fully polarized edges. While 
this indicates a mixture of ID and 3D features, area II can be 
identified as a quasi-lD phase because of the spatial pattern 
of the FFLO oscillations emerging at the trap center Area II 
is only found at interchain couplings t± < 0.3. Third, area 
III has a two-shell structure where the FFLO-type oscillations 
reside in the entire area of the partially polarized core, sur- 
rounded by fully polarized edges. Area III is found across 
the whole range of interchain couplings at intermediate-high 
polarizations below the transition to the normal phase. 

We find that the FFLO character evolves very differently 
with polarization in the quasi-lD and quasi-3D regimes. Fig- 
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FIG. 1. (Color online) Phase diagram of the 3D coupled-chain 
Hubbard model. The particle densities n-f4 and the order parameter 
A are calculated as a function of the interchain coupling and the 
polarization at zero temperature, (a) The oscillation amplitude Ao 
of the order parameter and the density difference — ni at the trap 
center. The phase at the trap center can be divided into three: the 
fully-paired superfluid (SF) (iif — ni, A 7^ 0), FFLO (oscillating 
n and A), and normal (A = 0) phases, (b) Phases I-III associated 
with the shell structures in the trap, explained in the text. 



ure|2]shows the shell structures with increasing polarization P 
at both sides of the 1D-3D crossover In the quasi- ID regime 
[Fig.[2|a)], the evolution is dominated by the expansion of the 
FFLO core. In contrast, in the quasi-3D regime [Fig.|2]^b)], the 
fully paired core shrinks with increasing P while the FFLO- 
type oscillations at the shoulders move toward the trap cen- 
ter. Near the crossover, these quasi- ID and quasi-3D features 
coexist. At the interchain coupling t± — 0.3, the oscilla- 
tions of the order parameter A become significant at both the 
trap center and edges, and the fully paired shoulders decrease 
from both sides as P increases. In addition, when P is in- 
creased to enter the area III, the oscillations of A along the 
trap show a distinct feature: At small t±, the amplitude of A 
becomes spatially uniform, despite the trap, at intermediate P 
[see Fig. |2|c)]. This can be compared to the zero derivative 
of the FFLO momentum as a function of the chemical po- 
tential.t^ The uniform oscillations occur only in the quasi- ID 
regime, raising the possibility that a peak signal of the FFLO 
momentum is more visible in this regime. 

Our phase diagram suggests that the optimum spot for ob- 
serving the FFLO state is extended over a significantly large 
area of the dimensional crossover In particular, the largest 
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FIG. 2. Evolution of the FFLO oscillations in the 1D-3D crossover. The profiles of the particle densities n-1-4 and the order parameter A 
along the chain sites i are presented for the two regimes of interchain couplings (a) t± — 0.2 (quasi-lD) and (b) t± — 0.5 (quasi-3D). In the 
quasi-lD regime, for small polarization P — 0.04, the FFLO-type oscillations at the center are surrounded by the fully paired shoulders. The 
region of oscillating A develops from the center, expands toward the edges as P increases, and then emerges over the whole area for P — 0.13. 
On the contrary, in the quasi-3D regime, the oscillations are initially at the partially polarized intermediate regions and spread toward the center 
as P increases. At finite interchain couplings, the far-edges of the trap are always polarized. In (c), the evolution of the oscillation envelope of 
A is shown with increasing P att± = 0.2 in the quasi-lD regime. Uniform oscillations occur along the trap at an intermediate polarization. 



range of polarizations associated with the FFLO-featured ar- 
eas (IIh-III) at the trap center is found around interchain cou- 
pling t± ^ 0.4. Since we still have two control parame- 
ters, the onsite interaction and the temperature, the phase di- 
agram should be further generahzed to discuss ultracold gas 
experiments which are conducted at finite temperatures and 
with a tunable atomic scattering length. We have examined 
also a stronger interaction U = —4 at small t±'s, includ- 
ing the ID limit (t± — 0) where the 3D scattering approach 
fails. At small interchain couplings, we have found that the 
main effect of a stronger interaction is to increase the criti- 
cal polarization of the transition to the normal phase while 
the change in the phase II area is relatively small. In the 3D 
limit, stronger interactions tend to keep the fully paired SF 
core even at larger polarizations, as observed in the previous 
experiments in the BEC regime,'^ which can make the 3D- 
FFLO area much smaller 

The stability of the FFLO state at finite temperature is a 
crucial question. We have examined the behavior at a finite 
temperature T = 0.05 near the dimensional crossover at the 
interchain coupling t± = 0.3. A low temperature algorithm is 
used!^ We find that finite temperature can stabilize the polar- 
ized superfluid (pSF) state at small polarizations while it de- 
stroys the quasi- ID FFLO character Figure [3] shows compar- 
isons between the density and order parameter profiles at two 
different temperatures T — 0.05 and T = 0. These profiles 
indicate that area II with the FFLO core melts at T = 0.05 into 
the pSF state, similar to the results in ID.^^ In the case of shell 
structures, the polarization from the FFLO area is easily redis- 
tributed in the trap at finite temperature to create a BCS-type 
order parameter with extra majority particles accommodated 
as thermal quasiparticles. In contrast, it turns out that area III 
with a trap- wide FFLO character at higher polarizations is not 
affected by the finite temperature examined. 

A comparison with the previous mean-field results in ID 
and 3D lattices shows the drastic effects of local quantum 
fluctuations: (1) Our results indicate that the large size of the 



FFLO area predicted by mean-field theory in 3D lattice^SE 
may have been overestimated. In the 3D limit, we find that 
the FFLO state is broken near local polarization 0.35 at a trap 
center that is approximately a quarter filled, which contrasts 
to the mean-field values 0.6 (|C/| = 5.14)** and 0.75 
i\U\ = 6)' at the same filling. This enhancement of the nor- 
mal state can be understood in light of the fact that including 
the particle-hole channel was shown to reduce pairing signif- 
icantly in lattices;^** DMFT includes full local quantum fluc- 
tuations causing such higher order effects. (2) In ID lattices, 
mean-field theorji^Sl predicted only a 3D-Iike shell structure 
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FIG. 3. Polarized superfluid phase at finite temperature. For in- 
terchain coupling t± = 0.3, the profiles at temperature T — 0.05 
are compared with those at T = for the particle densities n-1-4 and 
the order parameter A along the chain sites i. At polarization (a) 
P ~ 0.09, corresponding to phase II, the structure of the FFLO core 
surrounded by fully paired shoulders found at T = is completely 
changed into the polarized superfluid phase with a uniform order pa- 
rameter at T = 0.05. The order parameter oscillations at the edges 
are still observed at T = 0.05. In contrast, at a higher polarization 
(b) P ~ 0.14, corresponding to phase III, similar FFLO characteris- 
tics are identified at both temperatures T = 0.05 and T = 0. 
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with a fully paired center and polarized edges. In contrast, we 
find the reversal of the shell structure, similar to what has been 
predicted by continuum mean-field studyP^ (3) In addition, 
our characterization of the FFLO state does not find the well- 
separated domain wall between the sign changes of the or- 
der parameter predicted by the mean-field theories in latticed 
and in continuum!^ Our findings, the drastic decrease of crit- 
ical polarization, the shell structure reversal, and the absence 
of domain walls, emphasize the importance of local quantum 
fluctuations in lattices, regardless of the dimensionality. 

Our findings are directly applicable to future experiments. 
Recent experiments realized a weakly coupled 2D array of 
ID tubes with ultracold ^Li gases, observing the density pro- 
file characteristics of attractively interacting spin-polarized 
ID Fermi gases. This system can be further extended toward 
the quasi-lD regime by adjusting the optical lattice potentials, 
and if needed, discreteness along the tube direction can be re- 
alized by 3D lattices. Although it is nontrivial to experi- 
mentally reveal the oscillating order parameter, several meth- 
ods have been recently sugg ested to detect the FFLO state par- 
ticularly in ID systems, se^^^^and references therein. Our 
DMFT calculations have shown that the FFLO phase occupies 
a significant area of the phase diagram throughout the 1D-3D 



crossover. In addition, our calculations at finite temperature 
have found that the FFLO character is more stable at inter- 
mediate rather than small polarizations where the phase be- 
comes the polarized superfluid at low temperature. Our find- 
ings will help identify the presence of the FFLO state in ultra- 
cold atomic gases. In a general perspective, by using beyond- 
mean-field calculations, we have investigated the dimension- 
ality effects on the existence of the FFLO state in 1D-3D lat- 
tice systems. This may potentially help to understand also the 
puzzling repulsive Hubbard model counterparts such as the 
striped phase. 
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